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ABSTRACT 



Compact groups, with their high number densities, small velocity dispersions, and 
an interstellar medium that has not been fully processed, provide a local analog to 
conditions of galaxy interactions in the earlier universe. The frequent and prolonged 
gravitational encounters that occur in compact groups affect the evolution of the con- 
stituent galaxies in a myriad of ways, for example gas processing and star formation. 
Recently, a statistically significant "gap" has been discovered mid-infrared (MIR: 3.6- 
8fim) IRAC colorspace of compact group galaxies. This gap is not seen in field samples 
and is a new example of how the compact group environment may affect the evolution 
of member galaxies. In order to investigate the origin and nature of this gap, we have 
compiled a larger sample of 37 compact groups in addition to the original 12 groups 
studied bv lJohnson et al.1 (|2007l ) (yielding 174 individual galaxies with reliable MIR pho- 
tometry). We find that a statistically significant deficit of galaxies in this gap region of 
IRAC colorspace is persistant in the full sample, lending support to the hypothesis that 
the compact group environment inhibits moderate SSFRs. Using this expanded sam- 
ple, we have more fully characterized the distribution of galaxies in this colorspace and 
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quantified the low density region more fully with respect to MIR bluer and MIR redder 
colors. We note a curvature in the colorspace distribution, which is fully consistent 
with increasing dust temperature as the activity in a galaxy increases. This full sample 
of 49 compact groups allows us to subdivide the data according to physical properties 
of the groups. An analysis of these subsamples indicates that neither projected phys- 
ical diameter nor density show a trend in colorspace within the values represented by 
this sample. We hypothesize that the apparent lack of a trend is due to the relatively 
small range of properties in this sample, whose groups have already been pre-selected 
to be compact and dense. Thus, the relative influence of stochastic effects (such as the 
particular distribution and amount of star formation in individual galaxies) becomes 
dominant. We analyze spectral energy distributions of member galaxies as a function of 
their location in colorspace and find that galaxies in different regions of MIR colorspace 
contain dust with varying temperatures and/or PAH emission. 

Subject headings: galaxies: evolution — galaxies: interactions — galaxies: clusters - 
galaxies: statistics — infrared: galaxies 



INTRODUCTION 



The majority of galaxies reside in poor groups of galaxies ([Mulchaevl 120001 ) . With number 
densities similar to those seen in the centers of rich clusters, compact groups constitute a high- 
density sub-category of poor groups (jHicksonl Il982l ). This environment is thought to provide a 
local analog to conditions in the early universe when the average galaxy density was higher and 
interactions were more common and prolonged. 

Despite their similar densities, compact groups exhibit marked differences from the centers of 
clusters. This is especially true with regards to the interstellar medium (ISM), a crucial component 
in the process of galaxy evolution. For instance, the centers of clusters tend to have no neutral gas, 
while compact grou ps typically contain at least som e neutral gas, though they are known to be H I 
deficient as a class (jVerdes- Montenegro et al.ll200ll ). These H I deficiencies in compact groups are 
likely related to the ISM behavior in the outskirts of clusters, where t he degree of H I deficiency is 
correlated with the radial distance of the galaxy to the cluster's center (jGiovanelli &: Haynesl ll985). 
In addition, both compact groups and clust e rs show a correlation between H I deficiency and X- 
ray detectability ( Giovanelli k, Havnes 1985 : Verdes- Montenegro et alJboOll ). suggesting that the 
neutral gas undergoes a phase change via ionization and heating. 

A key difference between compact group s and galaxy cluster s is that galaxies in the cores of clus- 
ters typically have negligible star formation (jBalogh et al.lll998l : Ivon der Linden et al.ll2010h. while 
some compact groups can h ost moderate or even intense star formation (llglesias-Paramo &: Vflchez 
1999l : iGallagher et al.ll2010l ). However, IRAS observations of compact groups indicate normal levels 
of thermal infrared emission, suggesting that the star formation rates (SFRs) in these groups are not 
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genera lly enhanced (jAllam et alJll995l : IVerdes- Montenegro et al.lll998l ). In fact, Ide la Rosa et al 
(|2007l ) studied the stellar populations of elliptical galaxies in compact groups and find evidence of 
a mechanism that quenches star formation. This is again si milar to the outer regions of clusters, 
which exhibit suppression of SFRs out to three virial radii (jLewis et al.ll2002l ). Gas processing is 
also different in cluster co res and compact group s, where ram- pressure stripping doe s not appear to 
dominate gas processing (IRasmussen et al.ll2008l ). In addition. ICortese et al.1 ()2006l ) have studied a 
compact group falling into a cluster, and found that the galaxies have undergone preprocessing due 
to the compact group environment. Despite their high densities, it appears as though the compact 
group environment has more in common with the outer regions of clusters than with cluster cores. 



Recently, I Walker et al.l (|2010i ) examined the distribution of galaxies from 12 compact groups 
in IRAC (3.6-8.0 /im) colorspace and found evidence for a statistically significant gap between the 
location of galaxies with colors consistent with normal stellar populations and galaxies with colors 



indic ative of star formation activity (previously identified by iJohnson et al.l 120071 ; I Gallagher et al 



2008). In contrast, no gap is seen in the colorspace distribution of the core of the Coma cluster, a 
control sample of field galaxies, or a sample of interacting galaxies. While this is not surprising, the 
varying distributions in IRAC colorspace indicate that these environments have had a significant 
impact on the evolution of their galaxies. However, a similar (though less pronounced) gap was 
observed in the outer regions of the Coma cluster. This is evidence that the gap is unique to 
regions of enhanced galaxy density (relative to the field) in which the neutral gas has not been fully 
processed. 

The discovery of the gap in IRAC colorspace naturally leads to the question of the origin of 
the gap, which has been intractable due to the small size of the previous sa mple. Initial results 
comparing the IRAC gap with specific SFRs (SSFRs; iTzanavaris et al.ll2010l ) of the same sample 



of gr oups show that the gap in IRAC colorspace corresponds to a gap in SSFRs (jWalker et al 



2010), suggesting a shared underlying cause. In order to better investigate the region of colorspace 
identified as the gap, we need to increase the number of observed compact group galaxies, so 
as to understand the characteristics of any systems that may be found in this region as well as 
characterize any correlation with physical properties along the sequence in colorspace. In order 
to further populate IRAC colorspace, in th is paper we p resent the colors of galaxies for a large 
sample of Hickson Compac t Groups (HCG: lHicksonlll982l ) and Redshift Survey Compact Groups 
(RSCG: iBarton et al.lll996l ) at redshifts of z < 0.035 with a full suite of IRAC data available from 
the Spitzer Space Telescope archive, yielding 49 compact groups including 179 galaxies. 

In this paper, w e adopt the following terms: the "original sample" is the 12 HCGs studied in 
Johnson et al.1 (|2007j) and subsequent papers. The "new sample" consists of the 21 HCGs and 16 
RSCGs not previously studied by this group. The "full" or "expanded sample" is the combined 
original and new samples. 
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2. DATA 
2.1. Samples 

2.1.1. Compact Groups 

Because the goal of this project is to understand galaxy evolution in the dense environment of 
compact groups, we began with the HCG and RSCG catalogs. These catalogs were compiled using 
different selection criteria, thus the groups may have slightly different characteristics; however the 
RSCG criteria were chosen to c reate a catalog with properties similar to the HCG catalog. The 



HCG catalog was compiled by iHicksonl (|1982l ) using photometric plates from the Palomar Sky 
Survey. The criteria to be classified as a compact group were N > 4 (four members within 3 
magnitudes of the brightest galaxy), 6n > 3#g (no other galaxies located within a radius 3 times 
the radius containing the nuclei of the members), and iig < 26.0 (the surface brightness must be 



high). The RSCG catalog, compiled by iBarton et al.l (jl9961 ). utilized magnitude limited redshift 



surveys to locate compact groups and determine their members. The criteria to be a member of a 
compact group were to have a radial velocity difference of AV < 1000 kms -1 and galaxy separation 
of D < 50 kpc between neighbors. Galaxy groups meeting these criteria were included in the RSCG 
sample. There was no isolation criteria for the RSCGs, and as a result many are embedded in larger 
structures. There are several groups that appear in both the HCG and RSCG catalog; naturally 
these groups have only been included once in the full sample. The RSCGs include groups with 
three or more members. The original criteria for HCGs was four or more members, but some have 
been trimmed t o three because one foreground or background galaxy was found to be projected 



onto the group (IHickson et al.lll992l ). Including the RSCGs in the sample will allow future analysis 



of the effect of varying group properties and environments on the member galaxies. 

For this study, the HCG and RSCG catalogs were sub-selected to only include compact groups 
at a low enough redshift (z < 0.035) such that the polycyclic aromatic hydrocarbon (PAH) features 
do not shift out of their rest-frame bands, which would make an analysis of their distribution 
in colorspace infeasible. We searched the Spitzer archive for groups meeting this criterion, and 
included groups for which all four channels of IRAC data are available. This yielded 21 HCGs (in 
addition to the original sample of 12 HCGs) and 16 RSCGs, bringing the total number of compact 
groups in the full sample to 49 (see Table [1]) , and increasing the number of galaxies from 42 to 
179. Of these 179 galaxies, 5 are saturated in one or more IRAC bands (HCG 16d, HCG56d, HCG 
90d, HCG 92b, and RSCG 4b), and are not included in the fits or statistical analysis presented in 
this paper. Five of the 16 RSCGs (21, 44, 66, 67, and 68) are known to be embedded in larger 
structures. We have undertaken our analysis both with and without these groups, and in most 
cases, the results do not change. 
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2.1.2. Comparison Samples 



As in IWalker et al.l (|2010l ) , we compare our compact group sample with samples of galaxies 
i n other environments. These samples are o ur approximation to a "field" sample - LVL+SINGS 



(Dale et al 



2009 



20071 ). interacting galaxies (Smith et al.l 120071) . and t wo samples from the Coma 



cluster - the center and the infall region ([Jenkins et al.1 [20071) • See IWalker et al.1 (120101 ) for a 
discussion of these samples. 

In order to compare samples of similar characteristics, we apply a luminosity cut to our sam- 
ples, illustrated in Figure [TJ We selected the lowest luminosity bin with more than two compact 
group galaxies as the minimum luminosity for all of our samples. Thus our luminosity cut is 
log (L4.5 [erg/s/Hz]) = 27.5. Figure [2] shows the Hubble types for our compact group sample 
against LVL+SINGS. As can be seen in this plot, the LVL+SINGS sample is dominated by late- 
type galaxies (which is to be expected for the field), while the compact group sample shows a 
bimodal distribution of morphological types, with a large number of early-type galaxies. 



2.2. Observations 

As this study utilizes archival data, the observational setups are heterog eneous. In addition to 



the Spitzer data, we also utilize data from the Two Micron All Sky Survey (jSkrutskie et al.l l2006) 
to plot IR spectral energy distributions (SEDs) to shorter A. 



2.3. Photometry 



The IRAC images were first convolved to the 8 //m PSF in order t o ensure a common pixel 
scale and resolution, then photometry was performed using SURPHOT (IReines et al.ll2008l ). The 
program determines apertures by finding a specified contour level in a reference image, then mea- 
suring the flux within the same aperture for each image. The apertures were determined using a 
combined image of all four IRAC bands (weighted by A -1 ), using a contour level of 1 — 2a. Several 
background annuli were used, with inner radii ranging from 2 — 2.5 x the radius of the aperture and 
outer radii ranging from 2.5 — 3x the radius of the aperture, using both the mode and resistant 
mean. Uncertainties were calculated from the standard deviation of the fluxes measured using 
different backgrounds. Extended source aperture corrections were made using the algorithms of T. 
Jarret10, and average corrections were 4%, 2%, 9%, and 14 % at 3.6, 4.5, 5.8, and 8.0 fxm, respec- 
tively. There are slight differences from ljohnson etaL ( 2007 ). likely due to variations in background 
measurements, though the new data are still consistent with the previous measurements. 



: http: //web. ipac . caltech.edu/ staff /jarrett/ irac/calibration/ 
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Fig. 1. — Histograms of L4 5 for the compact group sample (striped) overlaid on the comparison 
samples (solid). The dotted vertical line indicates the luminosity cut. 
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Fig. 2. — Distribution of morphological types (following lHaynes fc Giovanellilll984l ) for the compact 
group galaxies (striped) and LVL+SINGS galaxies (solid). As expected for a field sample, the 
LVL+SINGS galaxies are dominated by late-types, while the CG sample has a bimodal distribution. 
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Table 1. CG Sample 



Group 


R.A. a 


Dec a 




AOR(s) 




(J2000) 


(J2000) 


(Mpc) 




H02 


00 h 31 m 30!0 


+08°25'52" 


136.6 


rll249152,rll252224 


H04 


00 h 34 m 15!9 


-21°26'48" 


108.9 


r21926400 


H07 


00 h 39 m 23!9 


+00°52'4l" 


56.4 


rll249408,rll252480 


H15 


02 h 07 m 3910 


+02° 08' 18" 


90.8 


r26875648,r26880512 


H16 


02 h 09 m 31!3 


-10°09'31" 


52.9 


rll249664,rll252736 


H19 


02 h 42 m 45!l 


-12°24'43" 


57.3 


rll249920 


H22 


03 h 03 m 31!3 


-15°40'32" 


36.1 


rll250176,rll252992 


H25 


03 h 20 m 43!7 


-01°03'07" 


84.5 


r26785904,r26880768 


H26 


03 h 21 m 54!2 


-13°38'45" 


126.9 


r23027712 


H31 


05 h 01 m 38!3 


-04°15'25" 


54.6 


r 11250432 


H33 


05 h 10 m 47 s 9 


+18°02'05" 


105.1 


r23027968 


H37 


09 h 13 m 35!6 


+30°00'5l" 


95.5 


r23028224 


H38 


09 h 27 m 38!8 


+12°16 / 51" 


124.3 


r23026432 


H40 


09 h 38 m 54!5 


-04° 51' 07" 


96.9 


r23025920 


H42 


10 h 00 m 21?8 


-19°38'57" 


53.3 


rll250688,rll253248 


H47 


10 h 25 m 48!4 


+13°43'54" 


135.9 


r23026688 


H48 


10 h 37 m 45!6 


-27°04'50" 


37.7 


rll250944,rll253504 


H54 


ll h 29 m 15!5 


+20°35'06" 


28.1 


r23027456 


H56 


ll h 32 m 3916 


+52°56'25" 


116.0 


r23026944 


H57 


ll h 37 m 50!5 


+21°59'06" 


131.7 


r23027200 


H59 


ll h 48 m 26!6 


+12°42'40" 


54.1 


rll251200 


H61 


12 h 12 m 24!9 


+29°ll'2l" 


52.1 


rll251456 


H62 


12 h 53 m 08!l 


-09°13'27" 


54.9 


rll251712 


H67 


13 h 49 m 03!5 


-07° 12' 20" 


109.9 


r26876160,r26881024 


H68 


13 h 53 m 40!9 


+40°19'07" 


40.7 


r 4483584,r27600128 
r27600384,r27600640 


H71 


14 h ll m 04!6 


+25°29'06" 


131.8 


r23028736 


H79 


15 h 59 m ll!9 


+20°45'3l" 


68.0 


r23028992 


H90 


22 h 02 ir W0 


-31°55'48" 


32.3 


rll251968,rll253760 


H91 


22 h 09 m 10M 


-27°47'45" 


99.1 


r26877184,r26881536 


H92 


22 h 35 m 57?5 


+33°57'36" 


88.6 


r6011392 


H96 


23 h 27 m 58!3 


+08°46'27" 


118.5 


rl2297728 
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3. COLORSPACE 
3.1. Full Sample 

The IRAC colors of the full sample are presented in Figure [3j which shows that an underdensity 
of galaxies is persistent, though it has become more akin to a canyon than a gap. As the sample 
is now comprised of 179 galaxies, there are more galaxies with colors that place them within the 
canyon, though there is still a dearth of galaxies in this region relative to the number of galaxies 
with MIR-blue and MIR-red colors. The plot in Figure [3] shows the data with a linear fit, however 
a new result seen in this full sample of 49 compact groups is the curvature of the distribution in 
colorspace, which will be discussed in ^6.41 To quantify this curvature, we fit a quadratic to the 
data, also shown in Figure Thus, we can undertake our analysis with respect to the line or the 
curve, both fit to the data. 

For the statistical analysis, we considered the MIR colors in two ways: rotated so that the line 
shown in Figure [3] became the axis, and unwrapped so that the curve shown in Figure [3] became the 
axis. For both methods, the sample was shifted to have a mean value of zero. We have redefined 
the gap as the canyon using both sets of colors, as illustrated in Figure [H In both plots, we define 
the canyon to be where the histogram is less than half its median value and Cmir is < 1.0. The 
bounds of the canyon are indicated by the vertical dotted lines. The blue side of the canyon changes 
depending on whether we define it using the rotated or unwrapped color - the canyon determined 
from the unwrapped data is slightly smaller. However, the results of our statistical analysis do not 
change. 

The plots in Figure [5] show the results of the Kolmogorov-Smirnov (KS) test comparing the 
rotated and unwrapped color distributions with a model of a uniform distribution, where D is the 
maximum difference between the sample and the model, and a gives the probability that the model 
matches the data. While we do not expect any particular sample to have a uniform distribution, 
using a model of a uniform distribution as a baseline facilitates comparisons between environments. 
The grey region indicates the canyon, which is clearly manifested in this plot as the flat portion 
of the cumulative distribution function (CDF). The canyon is less extended in this larger data set 
than the gap seen in the original sample, but it is still present, and a is small enough that we 
can conclusively state that the compact group sample is not drawn from a uniform distribution 
in colorspace. The two methods of unraveling the MIR colors yield the same result in the KS 
test. Further analysis uses the unwrapped colors, because the canyon galaxies as defined by the 
unwrapped sample are more cleanly defined in colorspace (shown in Figure [6]). Excluding the five 
RSCGs that are embedded in larger structures (21, 44, 66, 67, and 68) increases a by an order of 
magnitude, but does not change the result. 
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Table 1 — Continued 



Group 


R.A. a 


Dec a 




AOR(s) 




(J2000) 


(J2000) 


(Mpc) 




H97 


23 h 47 m 22!9 


-02°19'34" 


88.3 


r26877440,r26881792 


H100 


00 h 01 m 2O!8 


+13°07'57" 


71.9 


r26877696 


R04 


00 h 42 m 49!5 


-23°33'11" 


89.0 


r 12342784 


R06 


01 h 16 m 12!4 


+46°44'20" 


69.1 


r26227712 


R15 


01 h 52 m 4r?4 


+36°08'46" 


64.3 


rl8646272 


R17 


01 h 56 m 2r?8 


+05°38'37" 


75.5 


rl8030336 


R2lt 


03 h 19 m 36!6 


+41°33'39" 


68.7 


x-10483456 


R31 


09 h 17 m 26 f ;5 


+41°57'18" 


29.5 


r23113472,r23113728 


R32 


09 h 19 m 50!8 


+33°46'17" 


96.3 


r4438784 


R34 


09 h 43 m 12!6 


+ 31°54'44" 


25.7 


r4439552,rl6300032 


R38 


10 h 51 m 46!7 


+32°51'31" 


27.2 


x-16304640 


R42 


ll h 36 m 51!3 


+19°59'19" 


93.5 


r25498880 


R44t 


ll h 44 m 0O!6 


+19°56'44" 


93.6 


r3858688 


R64 


12 h 41 m 33 ! ;2 


+26°03'56" 


73.7 


rl0532608 


R66t 


12 h 43 m 17 sg 


+13°ll / 47" 


13.9 


r 4477440,rl8288896 


R67t 


12 h 59 m 32 f ;8 


+27°57'27" 


106.4 


r3859456 


R68t 


13 h 00 m 10?7 


+27°58'17" 


98.3 


r3859712,r3859968 


R86 


23 h 38 m 34 ! ?4 


+27°01'24" 


117.6 


rl0926592 



a Group RA, Dec and distance (corrected for Virgo+GA+Shapley) 
taken from NED. The cosmology used was H = 73 km s" 1 Mpc" 1 , 

^matter — 0.27, 3-11(1 ^vacuum — 0.73. 

^RSCGs known to be embedded in larger structures. 
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Fig. 3. — The colorspace distribution of the 179 compact group galaxies c omprising the full sa mple; 
the shaded box highlights the location of the gap in the original sample (jWalker et al.ll2010l ). The 
asterisks indicate galaxies that are mildly saturated in one or more IRAC channels, these 5 galaxies 
are not included in further analysis. We performed two different fits to the data, shown above. 
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Fig. 4. — Histograms of the left: rotated colorspace distribution and right: unwrapped colorspace 
distribution as described in M3.ll In both plots, the dashed line indicates the median, the horizontal 
dotted line is half that value. The vertical dotted lines indicate the boundaries of the canyon, defined 
to be where the distribution is less than half the median value with an MIR color < 1.0. In the 
rotated sample, 5% of the galaxies fall in the canyon, while for the unwrapped sample, 3% of 
galaxies fall in the canyon. 

Table 2. Results of Statistical Analysis 



Compact Groups vs 


N 


NcGs 


D 


a 




Consistent? 


Uniform Distribution 




174 


0.30 


2.0 x 10" 


13 


reject 


LVL+SINGS 


93 


174 


0.47 


3.8 x 10" 


12 


reject 


Interacting 


31 


78 


0.41 


1.6 x 10" 


-3 


reject 


Coma Center 


114 


173 


0.36 


6.0 x 10" 


-8 


reject 


Coma Infall 


47 


173 


0.18 


0.16 




do not reject 



Note. — Results of the KS test for the compact group distribution against 
a model and the comparison samples, unwrapped using the quadratic fit to 
the compact group distribution. D is the maximum deviation (possible 
range: 0-1) between the two CDFs, while a gives the probability that the 
two samples were drawn from the same parent distribution (which depends 
on the number of galaxies in the two samples, N, and Nqq s ). Based on a, 
we have determined whether the two samples are consistent. 
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Fig. 5. — KS test for the left: rotated colorspace distribution and right: unwrapped colorspace 
distribution of the full sample against a model of uniform distribution. The maximum deviation D 
of the sample from the model is indicated by the vertical line. This large value of D means a low 
probability that the compact groups are drawn from a uniform distribution, 2 x 10 _11 %. 




Fig. 6. — Left: The colorspace distribution of the compact group galaxies, with symbol color 
illustrating whether a galaxy is blueward (blue), in (green), or redward (red) of the gap. Right: 
KS test comparing the original sample of 12 HCGs (dashed line) with the new sample of 37 CGs 
(solid line). These two samples are very similar, with a fairly high probability of being drawn from 
the same parent distribution, 18%. 
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3.2. Subsamples 

To investigate how the compact group environment affects galaxy evolution, we investigated 
the colorspace distribution of various subsamples. If a clear trend with physical properties were 
found, it would provide a clear picture of how environment affects galaxy evolution in compact 
groups. 



3.2.1. Original vs New 

To insure that the properties of the expanded sample are consistent with the original sample, 
we compare it with the original sample. The CDFs of the new and original sample are shown on 
the right in Figure [6j which illustrates that the new sample is consistent with the original sample. 
The difference in the gap seen in the original sample and the canyon defined by the full sample 
(highlighted by the shaded region) is clearly illustrated by the different ranges of the flat portion 
of the CDF - the canyon is narrower and slightly less pronounced than the original gap. 



3.2.2. HCGs vs RSCGs 

To determine whether the different selection criteria of the HCG and RSCG catalogs find groups 
with different properties, we compare groups from the two catalogs. The colorspace distribution of 
the HCG galaxies and RSCG galaxies are shown in Figure El We see that the HCGs and RSCGs 
occupy colorspace very similarly. 



3.2.3. Separated by Physical Properties 
As compact group galaxies occup y colorspace differe ntly than galaxies in other environments 



(shown in Table [2] and §3.31 as well as IWalker et al.l l2010l ) , we expect to see a trend with physical 



properties of compact groups, especially the properties which differentiate them from other environ- 
ments. Thus we binned the groups by projected physical diameter and projected physical number 
density and examined colorspace as a function of these properties, as shown in Figure Both of 
these properties could be an indicator of how frequently or intensely interactions occur, which could 
plausibly affect the triggering of star formation and/or transformation from activity to quiesence. 
However, as these figures show, there does not seem to be any trend with these properties. Further 
investigation of the CDFs of the subsamples reveals that the shape of the subsample's CDF does 
not correlate with either of these properties. This will be discussed further in §6.21 
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Fig. 7. — Left: The colorspace distribution of left: the 124 HCG galaxies and right: the 50 RSCG 
galaxies. 




Fig. 8. — The colorspace distribution of the full sample, broken into subsamples based on left: 
projected physical diameter and right: projected physical density. Neither of these plots show a 
trend, indicating that these parameters do not play a key role in determining a galaxy's location 
in colorspace within the range represented by the compact group sample. 
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3.3. Comparison Samples 

The colorspace distributions of our comparison samples are shown in Figure [9l As this figure 
illustrates, the LVL+SINGS sample spans almost the exact same region of colorspace as the compact 
group sample, though there is no evidence for a canyon. The interacting sample does not occupy 
the blue region of colorspace, indicating that this sample does not contain any galaxies dominated 
by stellar emission. In contrast, galaxies in the center of the Coma cluster predominantly fall in the 
blue region of colorspace, with only a few galaxies whose colors indicate activity. Like the Coma 
center, the infall region of the Coma cluster shows a concentration of galaxies with normal stellar 
colors, but also contains galaxies with colors indicative of activity and exhibits an underdensity 
of points in the same region as the canyon in the compact group sample. We performed two- 
distribution KS tests comparing these samples with the compact group sample; the results are 
given in Table [2j 

It is important to consider the colorspace distributions of the samples in the context of the 
morphologie s of their galaxies, given in Figure [2] for the compact group and field samples. As 



discussed in lWalker et al.l (12010 ). untangling the effects of morphology on MIR color is non-trivial. 
To some extent, morphology and MIR color are expected to track each other; e.g., late-type spirals 
are typically starforming and therefore expected to have red MIR colors. Therefore, it could be the 
case that the different morphological types that dominate each of the comparison samples is driving 
the KS test results. For example, the interacting galaxy sample has no E/SO galaxies by selection, 
and the morphology-density relation means that primarily quiescent E/SOs are found in the Coma 
Core sample. Unlike these two samples, both the compact group and LVL+SINGS samples span 
the range of morphological types. While MIR color tracks morp hology quite well in the compact 



group sample, this is not the case for the LVL+SINGS galaxies (IWalker et al.ll2010l ). Thus, while 
the dearths seen in both colorspace and morphological distribution for the compact group sample 
may be caused by the same evolutionary process, the causal relationship is unclear. 



4. COLOR-MAGNITUDE DIAGRAMS 

Galaxy evolution is tied to galaxy size and buildup. If the compact group galaxies show a trend 
with mass or activity, this would be revealed in a color-magnitude diagram (CMD). Comparison of 
the CG CMD with LVL+SINGS and Coma Infall (shown in Figures QZ3 QT] and provide insight 
into galaxy properties. All three samples show a trend for MIR-red galaxies to be more luminous. 
This is not surprising, as we would expect that galaxies with activity would have more warm dust 
and PAH emission, thus increasing their Ls.o- One notable feature is the dearth of intermediate- 
luminosity, MIR-blue galaxies in the field sample, though this could be a selection effect, specifically 
that the SINGS sample was chosen to include "interesting" galaxies, thus mid-luminosity, MIR-blue 
galaxies may not have been included. 
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Fig. 9. — Colorspace distribution of the comparison samples (black plus signs) overlaid on the 
compact group sample (grey dots). 
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Fig. 10. — L3.6 (left) and Lg.o (right) color-magnitude diagrams for compact group galaxies. The 
plus signs represent galaxies below the luminosity cut (as discussed in §2.1.21 and Figure [Q), the 
asterisks indicate saturated galaxies. The shaded region indicates the canyon in IRAC colorspace. 
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Fig. 11. — L3.6 color-magnitude diagrams of galaxies in left: the field sample of LVL+SINGS 
and right: the infall region of the Coma cluster. The plus signs represent galaxies below the 
luminosity cut (as discussed in §2.1.21 and Figure [1]). The shaded region indicates the canyon in 
IRAC colorspace. 



5. SPECTRAL ENERGY DISTRIBUTIONS 



In order to investigate the relative contributions of amount and temperature of dust in compact 
group galaxies as a function of colorspace, we have binned SEDs of individual galaxies as a function 
of location in colorspace, shown on the left in Figure [131 Since we are interested in the ISM, we 
extracted the contribution from stellar light by scaling a 5000K blackbody to the K-band emission. 
If the dust properties (e.g. temperature, composition, excitation) are the same in galaxies across 
colorspace, and the variation in MIR color is simply due to different amounts of dust and PAH 
emission, we expect the slope of the stellar-subtracted SEDs to be the same regardless of location in 
colorspace. If the different MIR colors are in fact dominated by varying dust properties, we expect 
the slope of the stellar-subtracted SEDs to change as a function of MIR color. As can be seen on 
the right in Figure [13] the stellar-subtracted SEDs do have clearly different slopes as a function 
of MIR color, which means that the range of MIR colors we see in the CG sample is not simply 
due to varying amounts of dust between galaxies. An analysis of the LVL+SINGS SEDs yields a 
similar result. Thus, the trend in IRAC colorspace can be attributed to the properties of dust in 
the galaxies, and the canyon reflects a range of moderate dust properties. 



6. DISCUSSION 

We have compiled a full sample of 49 compact groups comprised of 33 HCGs and 16 RSCGs. 
We see that the deficit (previously identified as a gap, see SJTJ) between galaxies with colors consistent 
with normal stellar emission (MIR blue) and galaxies with colors indicative of activity (MIR red) 
is persistent in the larger sample. We investigate physical trends as a function of colorspace and 
look at the structure of the colorspace distribution. We also compare the compact group colorspace 
distribution with samples of galaxies in other environments. We consider the SEDs of the galaxies 
to assess whether dust quantity or temperature changes as a function of colorspace. 



6.1. The Canyon 



The gap seen in the original sample of 12 HCGs (IJohnson et al.ll2007l : IWalker et al.ll2010l ) needs 
to be redefined in terms of the full sample presented here. While there is still an underdensity 
between galaxies whose colors are consistent with normal stellar populations and galaxies whose 
colors are indicative of activity, it is no longer the clear-cut gap seen originally. Now with better 
statistics, it is sparsely populated and thus more of a canyon. As shown in Figure [H we quantified 
the location of the canyon via a histogram of the colorspace distribution. As the full sample has 
~ 4x more galaxies than the original sample, there is naturally a larger absolute number of galaxies 
in the canyon, though the relative number of galaxies are comparable (2.4% of galaxies from the 
original sample, 2.8% of galaxies from the expanded sample). The larger number of galaxies in this 
region will allow us to learn about properties of galaxies that exhibit these intermediate colors. 
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Fig. 12. — Lg.o color-magnitude diagrams of galaxies in left: the field sample of LVL+SINGS 
and right: the infall region of the Coma cluster. The plus signs represent galaxies below the 
luminosity cut (as discussed in §2.1 .21 and Figure [1]). The shaded region indicates the canyon in 
IRAC colorspace. 




Fig. 13. — Binned SEDs of compact group galaxies, normalized to K-band, as a function of location 
in MIR colorspace. The yellow SED is comprised of canyon galaxies. The plot on the left shows 
SEDs with contributions from both stars and dust. The plot on the right shows stellar-subtracted 
SEDs, allowing us to study the dust. 
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6.2. Impact of Subsamples 

The larger number of galaxies available in the full sample allows us to look for trends in 
colorspace. Comparing the HCGs with the RSCGs shows that they are consistent with being 
drawn from the same parent distribution. Breaking the full sample into subsamples based on 
projected physical diameter or density does not reveal any trends in colorspace. This is a puzzling 
result - we hypothesize that compact groups are quantitatively different from other samples in MIR 
colorspace due to our analysis of the compact group and comparison samples, yet we do not see any 
trends with these properties within the compact group sample. It appears that projected physical 
diameter and density are not the dominant properties in determining the colorspace location of 
compact group galaxies over the parameter ranges of this sample. It may simply be that the 
compact group environment is pre-selected to be compact and dense, and thus does not have a 
large enough variati on in group diameters or densities to reveal trends in colorspace. However, 



Johnson et al.l (|2007l ) found that the original sample exhibits a trend in colorspace with group H I 
richness, so this seems to be the crucial property in colorspace distribution. Regardless of a group's 
diameter or density, star formation relies on the presence of cold gas. 



6.3. Environment 



In agreement with previous work (j Walker et al.ll2010l ). the colorspace distribution of the com- 
pact group sample is statstically different from a control sample of "field" galaxies, a sample of 
interacting galaxies, and galaxies from the center of the Coma cluster. However, we cannot rule 
out the hypothesis that the compact group sample is drawn from the same parent distribution as 
the sample from the infall region of the Coma cluster. This supports the hypothesis that there is 
something special about this environment - one of high galaxy density where the neutral gas has 
not typically been fully preprocessed (meaning str i pped, ionized, expelled, or otherwise modified 
through a galaxy interaction process) . ICluver et al.1 (|201ll ) found that galaxies in or ne ar the canyon 



show anamolous H2 excitation, similar to that seen in the shock in Stephan's Quintet (ICluver et al 



2010, HCG92;). This could be related to a rapid end to star formation in compact group galaxies, 
leading to the canyon observed in IRAC colorspace. 



6.4. SEDs 

The full sample revealed a curvature in colorspace not apparent in the original sample that 
was sparsely populated. The changing shape of the dust-subtracted SEDs as a function of MIR 
color indicate that the colors are not dominated solely by varying amounts of dust within the 
galaxies. Rather, the dust temperature distribution and/or PAH contributions change as a function 
of MIR color. We are unable to determine dust temperatures for these galaxies because the IRAC 
bands contain significant PAH features so we cannot disentangle the PAH emission from the dust 
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temperature, and since the SEDs are not turning over, we do not cover the peak of the emission. We 
have looked at IRS spectra of compact group galaxies in different regions of IRAC colorspace, and 
found that the spectra of the bluest galaxies decrease with wavelength and show no PAH features. 
The PAH features then appear fairly quickly as you consider redder galaxies. They seem to be fully 
formed by the canyon region, and do not change significantly as you look redward. The primary 
change redward of the canyon seems to be the slope of the spectra (Walker et al. in prep). This 
suggests that galaxies undergoing activity have different dust compositions than quiescent galaxies. 

6.5. Conclusions 

We have determined that the distribution of compact group galaxies in MIR colorspace still 
shows an underdensity in the canyon region of colorspace. This distribution is most similar to that 
of the Coma infall region, implying a similarity in environment. The distribution of colors is caused 
by varying dust temperatures and PAH emission rather than varying amounts of dust. However, 
there are still many unexplained properties of this expanded sample. 
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